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We investigate the combined effects of weak disorder and a two-dimensional (2D) optical lattice 
on the collective excitations of a harmonically trapped Bose-Einstein condensate (BEC) at zero 
temperature. Accordingly, we generalize the hydrodynamic equations of superfluid for a weakly 
interacting Bose gas in a 2D optical lattice to include the effects of weak disorder. Our analytical 
results for the collective frequencies beyond the mean-field approximation reveal the peculiar role 
of disorder, interplaying with the 2D optical lattice and interatomic interaction, on elementary 
excitations along the 3D to ID dimensional crossover. In particular, consequences of disorder on 
the phonon propagation and surface modes are analyzed in detail. The experimental scenario is also 
proposed. 
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I. INTRODUCTION 

The measurement of the frequencies of collective exci- 
tations has emerged as a fundamental and precise tool 
to investigate the quantum many-body physics of atomic 
Bose-Einstein condensates (BECs) in unprecedent detail 
[J-Q. It not only provides an excellent confirmation of 
the mean-field predictions, but also represents a very ef- 
fective method for probing the beyond-mean-field effects 
[E|-0] , particularly in a weakly interacting Bose gas whose 
condensate fraction is close to 100%. Theoretically, such 
a gaseous BEC's system can be well represented by a sin- 
gle macroscopic wave function, consequently allowing for 
clear hydrodynamic formulations that provide analytic 
or semi-analytic results in describing dynamic behavior 
of a BEC's system. 

Along this line, Stringari [j| pioneered in applying the 
hydrodynamic equations for superfluid to describe the 
ground and excited states of a magnetically trapped di- 
lute interacting Bose gas at low temperature, yielding 
predictions in excellent agreement with the experiment 
[lj. The same scheme was later generalized by Kramer 
et al. 8] to include the significant effects of optical lat- 
tice, following the spectacular experimental realization 
of optical lattices [J Q . Recently, the dramatic effects 
caused by disorder on the quantum properties of a BEC 
have attracted intensive attention [9j- For example, it 
has been observed that even a tiny amount of disorder 
in the confining fields leads to a fractioning of quasi- ID 
condensates in waveguide structures on atom chips [loj . 
One particular striking effect of disorder arises from Ref. 
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11] that disorder is more active in reducing superfluidity 
than in depleting the condensate even at zero tempera- 
ture (T = 0) . Investigations on collective modes of a BEC 
in a random potential using hydrodynamic scheme have 
been carried out by Falco et al. in Ref. [l2j in lattice-free 
case. In view of the possibility to control both disorder 
and a BEC in an optical lattice almost at will [l3| , an im- 
portant direction consists in investigating the collective 
excitations of a trapped BEC in the combined presence 
of disorder and optical lattice. 

Another significant influence of disorder on Bose gases 
is that it enhances quantum fluctuations and the beyond- 
mean-field corrections thereof to the equations of state. 
Such corrections in the free space arising from inter- 
atomic interaction have been predicted by Lee, Huang 
and Yang (LHY) [T^ . Yet measuring these corrections 
remains challenging, as they are usually too small to 
be observed. Pitaevskii and Stringari [6| first pointed 
out observing these beyond-mean-field effects in the fre- 
quency shift of collective excitations. Later in Ref. [ll|, 
Orso et al. suggested that these frequency shift could 
be further magnified by introducing optical lattice to en- 
hance correlations. Furthermore, Ref. [l2j showed that 
disorder also caused a frequency shift in collective modes, 
with a magnitude comparable to interatomic interaction. 
It is therefore expected that the combined presence of 
disorder and optical lattice might considerably enhance 
the quantum fluctuations, ultimately rendering the ef- 
fects beyond mean field more observable in the experi- 
ment. On the other hand, disorder has been shown to 
cause a departure from Kohn's theorem for the dipole 
mode [lH, a signature to distinguish the effects beyond 
mean field due to disorder from that due to interatomic 
interaction. An especially appealing question, therefore, 
consists in exploring the beyond-mean-field effects arising 
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from the interplay between disorder and optical lattice on 
the collective excitations of a trapped BEC. 



II. HAMILTONIAN FOR A TRAPPED BEC IN 
THE PRESENCE OF WEAK DISORDER AND A 
2D OPTICAL LATTICE 



Hence, we are motivated to launch a comprehensive 
investigation on a trapped BEC, taking into account of 
the combined effects of disorder, optical lattice and in- 
teratomic interaction. Against this background, a 2D 
optical lattice is especially favorable, because it induces 
a characteristic 3D to ID dimensional crossover in the 
quantum fluctuations of a BEC 15], even in the presence 
of weak disorder [l6| . Loading a BEC into a 2D optical 
lattice, therefore, provides an effective tool to study a ID 
BEC and 3D BEC through asymptotic analysis [H], [l6| . 



In this paper, we focus on studying the collective ex- 
citations of a harmonically trapped BEC in the presence 
of weak disorder and a 2D optical lattice at T — 0. Ac- 
cordingly, we extend the hydrodynamic equations for su- 
pcrfluid in a 2D optical lattice to include the effect of 
weak disorder. Using the developed hydrodynamic equa- 
tions, we calculate the collective frequencies beyond the 
mean-field approximation. Our analytical results reveal 
the peculiar role of disorder, when interplay ing with a 2D 
optical lattice and interatomic interaction, on the excited 
states of a trapped BEC along the 3D to ID dimensional 
crossover. In particular, the effects of disorder on sound 
velocity and surface modes are discussed in detail. We 
also propose possible experimental conditions to realize 
our scenario. 



The outline of this paper is as follows. In Sec. II, we in- 
troduce the Hamiltonian for a magnetically trapped BEC 
in the presence of a 2D optical lattice and weak disorder. 
We first examine the case in the absence of disorder in 
Sec. Ill where we follow the scheme in Q for general- 
izing hydrodynamic equations in Ref. [5] to include the 
effects of a 2D optical lattice. Thereafter in Sec. IV, 
we account for the presence of random potential and ex- 
tend the results in Sec. Ill to include the effects of weak 
disorder. Using these equations, we analyze in Sec. V 
the combined effects of disorder and a 2D optical lattice 
on the collective excitations of a trapped BEC. In par- 
ticular, we calculate the beyond-mean-field corrections 
to the collective frequencies. Detailed discussions on the 
sound velocity and low-lying modes as well as quantum 
behaviors along the dimensional crossover are presented. 
Finally in Sec. VI, we summarize our results and propose 
the experimental conditions for realizing our scenario. 



The N-body Hamiltonian describing the Bose system 
at T = has the form Q 



H — fiN = J dr¥(r) 



h 2 v 2 

2 to 



- t*+V ex t(r) + V ran (r) 



|#t (r) |r (r) 



*(r), 



(1) 



where \l/(r) is the field operator for bosons with mass m, 
fi is the chemical potential, N — J (r)^ (r) is the 
number operator, and g — 4irh 2 a/m is the coupling con- 
stant with a being the s-wave scattering length in the free 
space. In Hamiltonian (P), V ran {v) and V ext (r) respec- 
tively represent the random potential and the trapping 
potential. 

The external trapping potential V ex t (r) in Hamiltonian 
([T]) is generated by a superposition of a 3D harmonic 
confinement Vh (r) of magnetic origin and a 2D optical 
lattice V opt (r) modulated along the x — y plane 



V ext (r) = -m (uj 2 x x 2 + uj 2 y y 2 + uj 2 z 2 ) 
+ sE R [sin 2 (q B x) + sm 2 (q B y)] 



(2) 



where oj x , uj v and uj z are the frequencies of the harmonic 
trap, s is a dimensionless factor labeled by the intensity of 
laser beam and En — h 2 q 2 B /2m is the recoil energy with 
Hqs being the Bragg momentum. The lattice period is 
fixed by d = Tt/qB- 

Disorder V ran (r) in Hamiltonian (JlJ is produced by 
the random potential associated with quenched impuri- 
ties EE E3 



V ran (r) 



N imp 

E 

i=l 



vflr-ril), 



(3) 



with v(r) describing the two-body interaction between 
bosons and impurities, being the randomly distributed 
positions of impurities and Ni mp counting the number of 
IV To obtain the concrete form of the pair potential v(r), 
we need to investigate the scattering problem between a 
boson with mass m and a quenched impurity with mass 
AI in the presence of a 2D optical lattice [I^. Here, we 
restrict ourself to the conditions of a dilute BEC system 
in the presence of a very small concentration of disorder. 
Thereby, the potential v(r) can be expressed by an effec- 
tive pseudo-potential v(r) — gi mp 5(r) [TlJ . For the sake 
of notation convenience, one can write <?; mp = 2irh 2 b/m 
HHEl, E3- Here, b effectively characterizes the lattice- 
modified strength of disorder, while the reduced mass for 
the boson-impurity pair coincides with m when the mass 
of impurity is taken to be infinite due to its quenched 
nature [17], i.e. liniM^oo mM/ (M + m) = to. 
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III. MACROSCOPIC DYNAMICS OF A BEC 
LOADED IN A 2D OPTICAL LATTICE 

Let us first consider the absence of disorder (V ran = 0) 
in Hamiltonian ([TJ) . In this section we shall follow the 
scheme in Ref. [8J and reproduce corresponding results 
for generalizing the hydrodynamic equations of superfhud 
Q to include the effects of a 2D optical lattice at T = 0. 
The validity condition for the presented hydrodynamic 
formulation requires relatively large s where the inter- 
well barriers are significantly larger than the chemical 
potential, while quantum tunneling is still sufficient to 
ensure full coherence [1, |1) . 

Under above conditions, a formulation for the macro- 
scopic dynamics can be accomplished by considering the 
long- wavelength limit of Hamiltonian ([lj and the action 
functional thereof for the model system. To this end, we 
start from the ansatz for the condensate order parameter 
within the tight-binding approximation 0, [l5| 



*(r) = Y, W(x-l x d)W(y-l y d)fi x , ly (z)e 



(4) 



Here, Si xt i is the phase of the (l x ,l y ) component of 
the order parameter, while W(x — l x d), W(y — l y d) and 
fl at l (z) are real functions. In addition, the W(x — l x d) 
and W(y — lyd) are further assumed to satisfy the pe- 
riodic conditions, i.e. W(x — l x d) = Wq(x) and W(y — 
lyd) = Wo(y) where both Wo (a;) and Wo(y) are localized 
at the origin. A transcription to the continuum limit 
is achieved by coarse-graining through the replacement 
Y^,i x i y — > (l/^ 2 ) / dxdy. One thereby introduces the 
smoothed phase S (Si(z) — > S(rj) and a coarse-grained 
density Q 



n(r) 



1 



ft 



(5) 



lyd. 



with x = l x d, y 

The long-wavelength limiting form of the Lagrangian 
functional can then be obtained by substituting ansatz 
(U) into Hamiltonian ([1]), applying coarse- graining and 
proceeding as Ref. 



L = (-ih J ®*^®dr + H-fj,N) 



dr 



dS 

h— - 2t [cos(dd x S) 



cos (ddyS)] 



^{d z Sf + V ho (v) + g -n{v) 



where 

t = - I dxW (x) 



-^ x + V opt (x) 



n(r), (6) 



W (x-d) (7) 



relates to the tunneling rate between adjacent walls with 
Vopt(x) = sERs'm 2 (qBx) being the external optical lat- 
tice applied along x direction, and g — 4irh 2 d/m is the 



lattice- renormalized coupling constant with a = C 2 a and 
C = d J^, 2 WQ(x)dx. In addition, Eq. © is derived 
based on the Thomas-Fermi (TF) approximation [5( in 
the large N limit where one can neglect the quantum 
pressure terms originating from the radial term in the 
kinetic energy ||. 

With Eq. ([6]), the equations of motion are derived by 
requiring that the action A = J Q Ldt be stationary under 
arbitrary variations in the macroscopic density n(r) and 
the phase S(r) with fixed endpoints, i.e. SA/Sn = and 
SA/SS = 0. The resulting equations are [l9[ 



d ± [n sin (dd ± S)} + —d z [nd z S] = 0, (8) 



dn 2td 
~di + ~h 

flc fc2 

h^r + Sfi-2t cos (dd±S) + — (d z S) 2 = 0, (9) 
at ' 2m 

where Sfj, = Vho + gn — fi is the change of the chem- 
ical potential with respect to its ground state value 
16]. In Eqs. © and two notations have 

been introduced for simplicity: d±[nsm(dd±S)] = 
d x [n sin (dd x S)] + d y [n sm(dd y S)] and cos (dd±S) = 
cos (dd x S) + cos (ddyS). Accordingly, one can define the 
superfluid velocity v in a 2D optical lattice as Q 



J x{y) 



m 
H 



d x (y)S, 



d z S, 



where the effective mass m* is introduced by 



m 2mtd 2 
m* 



n 2 



(10) 

(ii) 

(12) 



which accounts for the increased inertia of the system 
along the direction of optical lattice [20j, |2l| . The ground 
state density within the TF approximation is obtained 
by putting 5fi = that yields [f| Ht : 



n T F(r) = - 
9 



1 / 9 2 

V - 2 m ( u x x 



i i 



2 2> 
■U1„Z 



(13) 



Eq. (|T3|) shows a fattened cloud in the presence of a lat- 
tice due to increased repulsive interaction characterized 
by 9- 

In the limit of small oscillations, Eqs. ([8]) and ([9]) can 
be linearized by substitutions of Eqs. (fit)]) and (|TTj) and 
decomposing n(r) = riTf(r) + Sn(r) and v(r) = 6v(r). 
The resulting equations read, 



— + V • [n TF Sv] = 0, 

h— + Sfj, = o. 

at 



(14) 
(15) 



With Eqs. (|T0|) and pT|) . one can combine Eqs. (fT4|) 
and (|15[) and consequently obtain the equation for the 
propagation of a density disturbance in the model system 
at T = 0, 



771- 



d 2 Sn 
W 



0. 



(16) 
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where we have introduced the notion 11 



A. Hydrodynamic equations 



m d 



m d d 



m* dx ' V fn* dy ' dz 



(17) 



Eq. (|16p bears formal resemblance with the corre- 
sponding equation in Ref. [f| for small-amplitude density 
fluctuation in the free space. Hence by assuming that the 
effect of optical lattice is captured by the effective mass 
m* and the renormalized interaction g, one can study the 
quantum behavior of a BEC as if the space is homoge- 
neous in the absence of harmonic trap, despite the ap- 
plication of periodic potential. In particular, in the limit 
8t ^> /i where the system retains a 3D behavior, Ref. 
[HJ has experimentally verified the validity of the above 
mass renormalization theory. In the opposite ID regime, 
on the other hand, the use of hydrodynamic equations is 
justified by the superfluidity of the quasi- ID tubes [23[ 
and by the avoidance of the Mott-insulator phase. 



IV. MACROSCOPIC DYNAMICS OF A BEC IN 
THE PRESENCE OF WEAK DISORDER AND A 
2D OPTICAL LATTICE 

In this section, we take into account of V ran (r) ^ in 
Hamiltonian (p} and investigate the effects of weak dis- 
order, combined with a 2D optical lattice, on elementary 
excitations of a trapped BEC at T = 0. To this pur- 
pose, we generalize corresponding hydrodynamic scheme 
in Sec. Ill to include the effects of weak external ran- 
domness. 

The disorder in this paper is associated with quenched 
impurities which generate a random potential V ra n{ r ) 
given by Eq. ([3]). For simplicity, we further assume that 
the randomness is uniformly distributed with density 
nimp = Ni mp /V and Gaussian correlated [ljj ■ Hence, 
the two basic statistical properties of disorder, in the 
presence of a 2D optical lattice, are the average value 
(Vq) = (jimpnimp and the correlation function 



(VkY-u) = —g\ mv n mp , 



(18) 



where g% m p = 2itK 2 b/m is the lattice-renormalized boson- 
impurity interaction, ni mp is the impurity concentration, 
and Vk = (l/V) J e r V ran (r)dr is the Fourier transform 
of Vran( r )- Here, the notation (...) stands for the ensem- 
ble average over all possible disorder configurations [28[ . 
The effect of disorder in a lattice can be characterized by 
two important parameters: the concentration of disorder 
k — "n-imp/n and the ratio of effective interaction strength 



R = 



(19) 



The validity conditions for following hydrodynamic for- 
mulation in a random potential require the correlation 
length of disorder to be much smaller than the healing 
length of superfluid [12] , in addition to those required in 
Sec. III. 



We start from Hamiltonian ([T]) with V ran (r) given by 
Eq. ([3]) within the white noise approximation [17] . and 
follow similar procedures in Sec. III. The resulting hy- 
drodynamic equation for the macroscopic phase S reads 



ao fc2 
h— + S (i - It [cos (dd x S) + cos (ddyS)} + — (d z S) 2 = 0, 

(20) 

where Sfj, = V) 1G + gn — \i is the change of chemical poten- 
tial with respect to its ground state value. Compared to 
Eq. (|9]| in the absence of disorder (V ran = 0), Eq. (|20|) 
takes similar form, except that the associated chemical 
potential \x here includes the dramatic effects of disorder 
on the system's ground state [IJ HH HH ■ The presence 
of disorder, therefore, only affects the dynamics of phase 
fluctuation indirectly through the equation of state. We 
will determine fi in the following section. Once Eq. (|2"0")) 
is solved for S, the superfluid velocity v can be obtained 
from Eqs. (fT0|) and (fTTj) . In the limit of small-amplitude 
phase fluctuation S = SS, Eq. (I2TH) can be linearized as 



H ^T + ^ = 0. 
at 



(21) 



Here, the hydrodynamical variables are assumed to be 
ensemble averages over all possible realizations of exter- 
nal disorder. By this assumption, the validity of Eqs. 
(|20p and ([21]) should be restricted to the self-averaging 
regime where the wavelength of the hydrodynamic modes 
is much larger than the correlation length of disorder [l2j ■ 
For the equation of total density n(r), the correspond- 
ing conservation law is unaffected by the presence of weak 
disorder. However, the random potenti al g ives rise to a 
normal component in the fluid (TJ, [TH, Il6l I2H ] which is 
pinned by quenched impurities and is neither a fluid nor 
dynamic. Thus, when a small perturbation 5v is ap- 
plied, only superfluid component responds to the probe 
and contributes to the induced current j = n s 5v with n s 
being the superfluid density [2(1 [H| , whereas the normal 
component remains stationary. Consistent with this pic- 
ture, a density disturbance Sn(r) from the equilibrium 
density profile is related to the perturbation 6v through 
the equation of continuity 



dSn 



V • [n s 6v] = 0. 



(22) 



Here, Sv is determined from SS through Eqs. (fTU|) and 
(fTTj) . Hence Eq. ([22) can be combined with Eq. (|2TJ to 
yield the hydrodynamic equation 



d 2 5n 



0, 



(23) 



where the notation V is given by Eq. (fT7[) . 

Eq. (12"3"[) is one of the main results of this paper. This 
equation is characterized by the disorder-depleted super- 
fluid density n s , the renormalized interaction g and an 
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affective mass m* accounting for the presence of optical 
lattice. For vanishing disorder, the whole system partici- 
pates superflow and Eq. (|23| becomes the corresponding 
equation for superfluid in Refs. (8l. [l5j: whereas for van- 
ishing optical lattice, our result recovers the correspond- 
ing one in Ref. [T3| ; while for vanishing optical lattice as 
well as disorder, Eq. (|2"5|) passes to that in Ref. Q for su- 
perfluid in harmonic traps. We emphasize that Eq. (|23[) 
has no similarity with the corresponding phenomenologi- 
cal equations in Landau's two- fluid theory (24|. The dis- 
order induced normal component is not dynamic. It's 
effect on the dynamics of superfluid is indirectly exhib- 
ited through the superfluid density n s = n — n n and the 
equation of state. 



B. 



Equations of state beyond the mean-field 
approximation 



Eq. (|23|) needs to be supplemented with prescriptions 
for the chemical potential /x(r), the ground state density 
n(r) and the superfluid density n s (r). In what follows, 
we derive expressions for these quantities based on the 
local density approximation (LDA) 



1. Chemical potential 

For the density profile that varies on a macroscopic 
scale, one can invoke the LDA for the chemical potential 



H(r) = ^i[n{r)] + V ho (r), 



(24) 



where ui refers to the chemical potential calculated for 
a BEC's system with Hamiltonian ([TJ in the absence of 
harmonic trap (Vho — 0). 

The fii can be determined from ui = dE g /dN, with 
Eg representing the ground state energy of the model 
system without harmonic trap. To calculate E g micro- 
scopically, we adopt Bogoliubov's theory [l6[ and ex- 
pand the field operators in Hamiltonian (fTJ) in the form 

*( x ) = Ekake~^ 2 0fc x {x)(t> ky (y), where <j} kx (x)(/) kv (y) 
corresponds to the lowest Bloch band of the BEC sys- 
tem under consideration. In the tight-binding ap- 
proximation, one can write <f> kjs (x) = e llkx w(x — 
Id) with w(x) = exp\-x 2 /2a 2 ]/'K 1 / 4 o- 1 / 2 and d/a ~ 
7rs 1 / 4 exp(-l/4^/i) ja Ujj. The dominate presence of 
condensate is taken into account through the approx- 
imation fik ~ v 7 nVSkQ + Sdk. Retaining only the 
quadratic terms in the excitations Sa^ and Sa^ from 
the condensate, the truncated Hamiltonian ([I]) can be 
ultimately diagonalized by the Bogoliubov transforma- 
tion [11]: ak = u k c k - WkcL k - ^f~NV k («k - f^k) 2 /E k 



and 

2 
k 



V; 



MkCk - WkC_ k - 

- u k c-k - VNV-h (uk - Wk) z /Ek with 

= [(eP k + gn)/Ek-i\/2 and £ k = 

2gn). Here, Vk is the Fourier transform of 
h 2 kl/2m + 2t[2 — cos(k x d) — cos(k y d)] is 



u*4 
. - 1 



V4H 

V ran (r) and e£ - ,„ 
the energy dispersion of a noninteracting Bose gas in the 



presence of a 2D optical lattice. Consequently, one finds 
the ground state energy E g = gn 2 V/2 — X^o( £ k + 9 n ~ 

E w )/2 + N \g imp n imp - (g 2 mp n imp /V) Y,^o £ V E l ■ 

The resulting /i/ in Eq. (|24p . in a convenient form, is 
given by 



Hl(ri) = gn [1 + k in t(n) + k dls (n)] 



(25) 



which includes the first order correction to the result 
fi = gn in the mean-field approximation. Precisely, in- 
teratomic interaction and weak disorder respectively give 
rise to the beyond-mean-field corrections 



kint(n) = 



1 y/2mgn 
Airh nd 2 



df(x) 
dx 



;/(*) 



and 



b R yjmgn 
a AttH nd 2 



(26) 



,(27) 



L vm P 



jn. Here, the f(x) and 



with x = 2t/gn and k 
Q(x) as functions of variable x are respectively defined 
as uM 



fix) 



* d 2 k 2-Fi 



2y^ J-„ (2ir) 2 



(28) 



and 



Q(x) 



P I 1 1 1 xj 

7T /" dfcj; 12' 2' 1 ' (l+z+z sin 2 (fc !c /2)) 2 / 



27T 



^/xsin 2 (fc^/2) + I 



(29) 

In Eq. d2HJ) , l(k) = 2 - cos(fc x ) - cos(fc y ). The function 
2Fi(a,b,c,d) in Eqs. (f28f and (|29|) is the hypergeomet- 
ric function [3ll ] and the integration over the transverse 
quasimomenta is restricted to the first Brillouin zone, i.e. 
\k x \, \k y \ < tt in Eq. (J28]) and \k x \ < tt in Eq. (|29|. 

The parameter x = 2t/gn controls the dimensional 
crossover for a uniform gas (Vtr = 0). In the limit x 0, 
corresponding to 8t -C ^; , the system undergoes a dimen- 
sional crossover to a ID regime where the f(x) saturates 
to the value 4\/2/3. In this limit, we can neglect the 
Bloch dispersion and Eq. (|25|) approaches asymptoti- 
cally to the chemical potential of a ID Bose gas in the 
presence of weak disorder 



Hi(n) = gn{ 1 



I \Jmgn 
nh nd 2 



b R y/mgn 
a Aith nd 2 



dQjx) 
dx 



\Q{x) 



'•(30) 



Eq. (|3T))) effectively generalizes the Lieb-Linger solution 
of the ID model expanded in the weak coupling regime 
to include the effects of weak disorder [l6j . 

In the opposite limit x ^> 1, corresponding to 8t^ui, 
the system retains an anisotropic 3D behavior. In this 
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regime, the functions (|28|) and (|29|) respectively reach 
the asymptotic law f(x) ~ lA3/y/x — 16v2/15to and 
Q(x) ~ —l/y/nx. Eq. (|25]l thus takes the asymptotic 
form 



1 



32 to* 
3i/7r m 



(31) 



In Eq. ([3"Tjl. the («i>)/5 and a/a cr = -(0.24a/ 'd)y/rrF/m 
[l5| provide a further renormalization of the scattering 
length due to the random potential and optical lattice, 

i.e. g e f f = g (l + nb/a + a/a cr ^j ; whereas, the remaining 

two terms proportional to the renormalized gas parame- 
ter generalize the LHY correction [l4[ to include effects 
of optical lattice and disorder. Eq. (|31l) therefore shows 
that the effects beyond mean field are much amplified via 
the introduction of periodic potential and external ran- 
domness, due to the increased inertia of mass m* along 
the direction of lattice and the interaction of bosons with 
impurities characterized by R. For vanishing disorder 
k = 0, Eq. (|3"T|) is reduced to the corresponding result 
in Ref. [1 5j j ; whereas for vanishing optical lattice s = 0, 
our result recovers exactly the corresponding one in Ref. 



17]. 



2. Ground state density profile 

Using Eq. (|25|) . one obtains by iteration the LDA 
ground state density profile n(r) [5[ 

ra(r) = n TF - n TF k int (n TF ) - n TF k dis (n TF ), (32) 

where 

n TF (v) = n(0) - fr(<^ 2 + u>y + u? z 7?) (33) 

is the mean-field value of density within the TF approxi- 
mation. In Eq. (13"3"|) . n(0) is the lattice-renormalized TF 
density evaluated at the center of harmonic trap deter- 
mined from 

n(0) = g(^ 2 +^y 2 +^ 2 ), (34) 

where X, Y, Z denote the size of the system. The second 
and third terms in Eq. (|32l) show that the quantum fluc- 
tuations due to interatomic interaction and disorder re- 
duce the LDA total density, thereby fattening the cloud. 

In the asymptotic 3D regime, k% n t in Eq. (|26[) and kdis 
in Eq. (|2"T|) respectively approach 



and 



, , a 32 to* 3 

hm (r) = 1- ^— f= — (Tia) 2 , 

a cr oy/ir to 



b ~ to* 

kdisij) = K- + 5R (na) 

a m 



(35) 



(36) 



It thus follows from Eq. (|32|) that the total density in 
the 3D limit takes the form 



n(r) = n FF 



b 

■K— - 

a 



( 32 



+ 5R 



TO 3 

{an TF ) 2 

TO 



(37) 

Eq. (|3"T)l shows that the effects beyond mean field due to 
both optical lattice and disorder have two consequences 
on the total density: first of all the TF density is further 
renormalized arising from the modified effective scat- 
tering length n FF (l ~ nb/a — a/dcr^j ~ griT F /g e f / = 

n e ff] while the second point concerns the generalized 
LHY correction in the presence of disorder and optical 
lattice as exhibited by the last term in Eq. (|3"T)l . For van- 
ishing disorder and optical lattice, Eq. (|37"f is reduced to 
the corresponding one in Ref. @. 

It's important to mention here that the density n(r) 
entering the hydrodynamic equations is the total density 
and shall not be confused with the condensate density 
n c (r) [6j, particularly in cases where the effects beyond 
mean field are taken into account. In the asymptotic 3D 
regime, by extending the result of quantum depletion in 
16] for an effectively uniform gas to that with harmonic 
trap using LDA, and together with Eq. (|3T|) for the total 
density, one obtains the condensate density of the model 
system 

n c (r) = n TF [ 1 - K- - J^- ) 
\ a a cr 



40 . v^F x - 
3^ + (5 + ~ )R 



(&n TF )*. (38) 



Comparison of Eqs. ([3T]) and (f3"5]) shows that the effects 
beyond mean field reduce more condensate density than 
the total density. The n(r) and n c (r) only equal in the 
limit where quantum depletions are negligible. In par- 
ticular, Eq. ([3"5]) demonstrates that the disorder induced 
quantum depletion 

, % b , \/tt „ - TO* , .3 

n R (r) = K-n TF + (5 + )R (an TF ) 2 (39) 

a 2 m 

is considerably amplified at each point of the space, com- 
pared to the corresponding result in Ref. [ll[ , when 
beyond-mean-field effects are included. 



3. Superfluid density profile 

The superfluid density is determined from n s (r) = 
n(r) — n n {r). Here, n(r) is the total density determined 
from Eq. (|3"2"1) . while n ra (r) is the normal fluid density. 
For a uniform gas in the absence of harmonic trap, the 
normal density n n can be determined from the long- 
wavelength limit of the static transverse current-current 
response function xr(q) at T = [TTI. l27t IM l33lj 



n n = limxr(q)) 



(40) 
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where xt(<1) can be derived within Bogliubov's theoreti- 
cal framework [ll], E3 • The result for n n of a BEC in the 
presence of weak disorder and a 2 D op tical lattice has 
been derived in our previous work [1 6j j within the con- 
text of renormalized mass theory. Thus by extending the 
corresponding result in Ref. [16] to that with harmonic 
trap using LDA, we obtain 



n„{r) 



R 

AHcP 



21: 



y / 2mgn T Fl{- ) 



griTF 



(41) 



Here, the «tf is the position-dependent TF density given 
by Eq. ([55)1 and the function I(x) with x — 2t/griTF is 
defined as 



I(x)- 



d 2 k 



1 



(27T) 



fx + y/xj(k) + 2 

(42) 

where 7(k) = 2 — cos^.) — cos(k y ). Here, the parame- 
ter x — 2t/griTF controls the dimensional crossover for 
nonuniform gases. In the limit of 3D regime correspond- 
ing to i > 1, I(x) asymptotically approaches the func- 
tion \/2/(6irx). Accordingly, the LDA normal fluid den- 
sity (|41l) takes the asymptotic form 



2y / 7r ~m* , „ 3V 

n n (r) = — — R n TF ln T Fa ) 

6 m ' 



(43) 



However, with the disorder induced quantum depletion 
given by Eq. ()39j) . the ratio n n /nR at each point of the 
space departures from the factor 4/3 pointed out in Ref. 
[ri| . as a result of the effects beyond mean field. 



V. COLLECTIVE EXCITATIONS 

Having at hand all the quantities needed to solve Eq. 
(|23|) . we can now proceed to investigate the low-energy 



collective excitations of the BEC's system under consid- 
eration. For this purpose, it's more advantageous to write 
n s = n — n n and Sfi — (dSfii/driTF)8n in Eq. (|23| with 
[i\ given by Eq. ([23)) . After some elaborated algebra, we 
ultimately obtain 



mu! 2 Sn + V • 



riTF^gSn 
( - 2 dkdis 



= -V 
5n | + V 



~ 2 dki n t 

9 n TF~n dn 

on TF 
n n VgSn 



(44) 



Eq. (|44[) provides the appropriate generation of the ze- 
roth order hydrodynamic equation in a 2D optical lattice 



/8n 



riTF^gfin 



= 0, 



(45) 



where the frequency oj can be solved for low-energy col- 
lective excitations within the mean-field approximation 
[1, Hj]. After solutions to Eq. (|4"5|) are found, one can 
then solve Eq. ([4"4"]l by treating all the terms on the right 
side as perturbation. Two types of such beyond-mean- 
field corrections can be identified: the first is associated 
with increased bulk compressibility x _1 = d/i/dn due to 
the enhanced quantum fluctuations caused by the com- 
bined effects of interatomic interaction, disorder and op- 
tical lattice [Hj]; the second type, by contrast, is unique of 
disorder which depletes superfluid fraction even at T = 0. 

Consequently, we follow Ref. @ and obtain the ana- 
lytical expression for the fractional shift in the frequency 
of low-lying collective modes 



5u> 



J dr\7 2 8n* 



TF 1 



dki. 



J dr V 2 Sn* (? 



dk di 



TF dr, 



n n VSn* ■ V5r 



2mu! 2 



J dr8n*5n 



J drSn*5n 



(46) 



where lo is the frequency of collective modes obtained 
from Eq. (|45[) within the mean-field approximation Q. 



Eq. (|46|) represents a major result of this paper. In 
what follows, we use this equation to derive some an- 
alytical results for phonon propagation and low energy 
collective modes. 



A. Phonon propagation 

An asymptotic analysis of Eq. (|4"6"]l allows us to ob- 
tain the sound velocity of the model system in the 3D 
limit. Particularly, when the harmonic trap is absent, 
the model system is effectively uniform within the con- 
text of mass renormalization theory [1, [l^- As a result, 
the titf, dkint/driTF, dkdis/driTF and n n in Eq. (|46|) are 
position-independent and can be taken out of the inte- 
gral. Together with Eqs. (1551) . (f5Bj) and (|4"3")l . one obtains 
the asymptotic law for the frequency shift in Eq. (|46[) in 
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the 3D regime 



Slu 



3 jl 

)(a + P + 7) u tf 



jdr (y 2 5n*^ 5r 



4mu> 2 



J drSn*Sn 



(47) 



Here, a = (32m* /Sy/jrm)a 3/2 , /3 = (5m* /m)Ra 3 / 2 and 
7 = — (4-^/77771* /3m) fla 3 / 2 . Meantime, solutions of Eq. 
(l4"5"j) for this effectively uniform Bose gas have the form 
5n ~ e lq r with a phonon dispersion u) — cq. Accordingly, 
we obtain the shift in the sound velocity 



5c 


* 

m 


' 8 
















c 


m 









(riTF<i 



s3\ 



(48) 



with respect to the lattice-modified Bogoliubov value 
c = y/gn/m*. The first term in Eq. (|48[) is a conse- 
quence of interatomic inter action, which generalizes the 
Beliaev result of 5c/ c — 8^na 3 /ir [HI for the sound ve- 
locity to include the effect of optical lattice. On the other 
hand, the remaining two terms (15/4— y / n)R(riTFa 3 ) 1 ^ 2 
are proportional to R, representing the influence of weak 
disorder on phonon propagation. As is shown, disorder 
affects the sound velocity in two opposite ways: the first 
is related to quantum fluctuations that tend to increase 
sound speed, whereas the second is related to the de- 
pletion of superfluid fraction at T = which tends to 
suppress the phonon propagation. The former effect is 
also shared by interatomic interaction, while the latter, 
by contrast, is a unique consequence of disorder. In ad- 
dition, all three terms in Eq. (|48jl demonstrate the same 
dependence on the gas parameter titf® 3 as the ground 
state chemical potential in Eq. (f3~Tj). Eq. (|4"5|) therefore 
represents generalized LHY correction to the sound ve- 
locity in the presence of optical lattice and weak disorder. 

It's worth mentioning here that one can also employ 
the compressibility-based definition for sound velocity 
[20L I2H I2H l35l - 437| to investigate the beyond-mean-field- 
effect induced shift. In our case, the normal compo- 
nent of fluid is pinned by the quenched disorder and 
doesn't directly participate in the propagation of the 
density disturbance. The corresponding definition for 
sound velocity typical of a superfluid in an optical lat- 
tice is therefore m*c 2 = n s d/j,/dn (26|. Here, the effect 
of disorder only enters indirectly through the compress- 
ibility x — dfi/dn and the expression for n s . To carry 
out the calculation beyond mean-field approximation, we 
consider the 3D asymptotic regime where the chemical 
potential approaches Eq. (|3"Tj) . The resulting sound ve- 
locity c s — c + 5 c is 




15 
T 



R 



-(riTFCi 



(49) 

The corresponding fractional shift is 5c/ c = [8/^tt + 
(15/4- y/^/3)R](m* /m)(n T Fa 3 ) 1/2 . This result, how- 
ever, is only consistent with Eq. (|48| for varnish- 
ing disorder (R = 0). In this case, both results re- 
cover corresponding ones in Refs. @ in free space and 



[15| in the presence of optical lattice. Nonetheless, in 
the presence of disorder, the two differ by an amount 
(2 v / 7rm*/3m)^(nTFa 3 ) 1 ^ 2 . This inconsistency in the 
presence of disorder may be attributed to the pathology 
of perturbation theory and its failure (at finite order in 
an expansion in powers of titfcl ) to deal with the to- 
tal density, condensate depletion and normal density in 
a completely consistent way. 



B. Collective modes 

Eq. (|46| allows some direct analysis of respective roles 
played by interatomic interaction and disorder on the 
elementary excitations of a trapped BEC. The effects be- 
yond mean field on low-lying collective modes have al- 
ready been analyzed in free space @, in the presence of 
optical lattic e and in the presence of disorder (but 
lattice free) [12[, respectively. The following analysis 
presents a modest generalization of these previous results 
to the combined presence of optical lattice and weak dis- 
order. 

According to Eq. (|4^|) . interatomic interactions will 
not affect the frequencies of the so-called surface modes 
that satisfy V 2 5n* — 0, despite the presence of optical 
lattice and disorder. In particular, the dipole mode corre- 
sponds to the center-of-mass motion of the superfluid and 
will not be influenced by the interaction, in accordance 
with the generalized Kohn's theorem 29]. Whereas, dis- 
order will have an immediate effect on surface modes, 
shifting corresponding mean-field frequencies in Ref. Q 
by an amount 
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J dvn n S75n* ■ V5r 



2mu 2 



J drSn*5n 



(50) 



In particular, the frequency shift in the dipole mode 
caused by disorder, as has been pointed out by Ref. [ijj], 
presents a departure from Kohn's theorem. Here, the 
result of Eq. ((511)) for dipole mode exhibits a magnified 
departure via the introduction of optical lattice. Such 
departure immediately distinguishes the beyond-mean- 
field effect due to disorder from that due to interatomic 
interaction. By contrast, the compressional modes will 
be affected by both disorder and interatomic interaction. 

To explicitly observe the effects of disorder, inter- 
atomic interaction and optical lattice on low-lying col- 
lective modes, we consider a deformed harmonic trap in 
the form 



2„2 



(51) 



where rj_ = (x^ + y 2 ) 1 / 2 is the radial coordinate and ui±_ is 
the radial trap frequency. The anisotropy is characterized 
by the parameter A = uj z /uj^_. 

For such anisotropic trap, the zeroth order hydrody- 
namic Eq. (14"5"T) admits solutions 5n(r) ~ r l Yi m (0 , 4>) for 
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m = ±1 and m = ±(Z — 1) [|[ that describes surface exci- 
tations, with f — \]m* /m(x 2 + y 2 ) 1 ^ 2 being the lattice- 
renormalized radial coordinate. The corresponding dis- 
persion laws are given by uj(m — ±1) — ylmjniFoJ±_ and 

oj(m = ±(l - 1)) = [(I - l)m/m*u\ + lu 2 ] 1/2 . The pe- 
riodic potential therefore effective ly mod ifies the radial 
trap frequency uj± by a factor \J m/m* , whereas leav- 
ing uj z unaffected. The dipole modes and quadrupole 
modes are of particular experimental relevance. The low- 
est dipole mode (I = l,m = 0) describes the oscillation 
of the center of mass along the z-direction 6n(r) ~ z, 
with the zeroth order dispersion given by w — uj z which 
is unaffected by the lattice. From Eq. ([50)) , the result for 
the disorder induced frequency shift in the dipole mode 



Suj_ _ 15 R 
~ ~ ~16~frcP 



1 2mg 



. p 2 \J\ - p 2 I . 
"(0) Jo V 1 - P 



I) 



dp, 
(52) 

where the variable y = 2t/gn(0) with n(0) given by Eq. 
(|34p and I(x) is a function of variable x — y/(l — p 2 ) 
given by Eq. (|42|) . Eq. (|52[) shows that the frequency 
shift is directly proportional to the strength of disorder 
characterized by R. For varnishing disorder (R = 0), 
there is no frequency shift in the dipole modes, in con- 
formity with Kohn's theorem (29|. For the quadrupole 
mode 5n(r) ~ f 2 sin 2 9e 2l< ^ with I = 2, m = 2, the mean- 
field dispersion is given by u) — y /, 2m/m*uj± . The corre- 
sponding frequency shift is found to be 



Soj _ 21 R 
~ ~ ~~8hcP 



1 2mg 



1 



r 



dp, 



which also varnishes in the absence of disorder, consistent 
with our previous analysis. 

Both Eqs. (|52l) and (|53p exhibit a dependence on 
the parameter y = 2t/gn(0) that characterizes the in- 
terplay between the kinetic energy in a lattice and the 
mean-field atom-atom interaction of a trapped gas. In 
the asymptotic 3D regime where y 3> 1, the func- 
tion I(x) with x = y/(l - p 2 ) in Eq. ([S3) de- 
cays as I(x) ~ s/2/6ttx. Accordingly, one finds 
the frequency shift in the dipole mode as Slo/lj — 
— (57r 3 / 2 m*/64TO)i?'\/a 3 n(0) and in the quadrupole mode 
as Su/u) = -(2l7r 3 / 2 m* /256m)R v / d 3 n(0). Both results 
exhibit dependence on the gas parameter and disorder 
strength in a similar way as the normal density does in 
Eq. (j43|). which indicates that the peculiar role of disor- 
der on the surface modes is closely related to the presence 
of normal component due to disorder. 

To observe the effects beyond mean field on the com- 
pressional mode, we further assume an effective disc- 
shaped trap ui z 3> \J m/m*ui± for simplicity. In this case, 
the zeroth order dispersion for the lowest compressional 
mode (n = 1, 1 = 0) is given by uj — y3u) z , correspond- 
ing to a density fluctuation 6n(r) ~ z 2 — 2p/3mu> 2 . Note 
that this compressional mode is only determined by the 



trap frequency ui z , thereby unaffected by the lattice ap- 
plied along the x — y plane. It follows from Eq. (|4T>)) that 
the frequency shift in the mode (n = 1, 1 = 0) is 



5u) 



1 2mg 



21 1 



512 hd 2 V n(Q) 



K(y)+R[(Z(y)-T(y)] . (54) 



Here. 



y = 2t/~gn(0), and K{y) = (l6/n)y^ 2 p 2 (p 2 - 
l)D{x)dp with x = y/(l - p 2 ) and D{x) = x 3 / 2 f"(x) - 
3f(x)/4^/x. Using the same definitions for vari- 
ables x and y, the functions Z(y) and T(y) are 
defined as Z(y) = {^/^y 1 ' 2 J* p 2 (p 2 - l)G{x)dp 
where G{x) = y 3 / 2 Q"{x) - 3Q(x)/Ay/E, and T(y) = 
32 $1 pV 1 - p 2 I{x)dp where I(x) follows Eq. (g2J|. Eq. 
(1541) confirms our previous statement that the compres- 
sional mode is affected by both repulsive atom-atom in- 
teraction and disorder. Specifically, the first term in Eq. 
(|54l) represents the contribution from interatomic interac- 
tion, whereas the latter two are consequences of disorder. 
Moreover, Eq. (l54l shows that disorder has two opposite 
effects on compressional excitations: the induced quan- 
tum fluctuation tends to increase excitation frequencies, 
whereas the depiction of superfluid fraction tends to de- 
crease the excitation energy. This latter effect exhibited 
by the term proportional to T(y), in particular, is unique 
of the effect of disorder on compressional modes. For 
vanishing disorder R — 0, Eq. (|54j) recovers the corre- 
sponding result in Ref. [15] . 

The parameter y = 2t/gn(0) controls the dimen- 
sional crossover of the frequency shift in Eq. (|54|) . In 
the asymptotic 3D regime corresponding to y ^> 1, 
functions f(x), Q(x) and I(x) with x — y/(l — p 2 ) re- 
spectively decay as f(x) ~ lA3/y/x — I&\f2/15-Kx, 
Q(x) ~ —I/^/ttx and I(x) = ^/2/Qnx. Hence 
Eq. ([Ml) asymptotically approaches Soj/uj = 

35V^rn*/128m [l - 3/5.R (250F/64 - vr/8)j y/dh^). 

This result provides a further generalization by disorder 
of the frequency shift caused by the lattice-generalized 
LHY correction in Rcf. Moreover, this result shows 
that the disorder-induced frequency shift is of the same 
order of magnitude, but with an opposite sign, compared 
to that due to atom-atom interaction. 



VI. POSSIBLE EXPERIMENTAL SCENARIOS 
AND CONCLUSION 

Central to test the validity of the hydrodynamic pic- 
ture proposed in this paper are the experimental abilities 
to measure the superfluid density and condensate density 
respectively. For a gaseous BEC, the condensate fraction 
is readily measured through mapping of occupation num- 
bers in the momentum space to real space by expansion 
imaging [lj. Whereas, to determine the superfluid frac- 
tion, one can use a method proposed in Ref. [38j . There, 
a BEC is beard to an optically-induced vector potential, 
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which simulates a uniform rotation of BEC. As a result, 
the normal fluid picks up non-zero angular momentum, 
whereas the superfluid acquires no angular momentum. 
Thus the spectroscopy can be used to measure the net 
change in the angular momentum of the fluid, and there- 
fore the superfluid fraction. 

Another difficulty that may arise is how to experimen- 
tally measure the beyond-mean-field correction to the 
equation of state due to the interatomic interaction and 
disorder. The concern is as follows [3^ : for a harmonic 
trap, the quantum depletion cannot be observed during 
the ballistic expansion in the typical Thomas-regime. Be- 
cause the mean-field energy is much larger than the trap 
frequency, the cloud remains locally adiabatic during the 
expansion. The condensate at high density transforms 
adiabatically into a condensate at low density with di- 
minishing quantum depletion. The above concern is true 
for a harmonically trapped BEC. However, it can be ruled 
out by introducing an optical lattice. In general, the con- 
finement frequency at each lattice site far exceeds the 
interaction energy, and the time-of-flight images are es- 
sentially a snapshot of the momentum distribution at the 
time of the lattice switch-off, thus allowing for a direct 
observation of the quantum depletion. 

Upon overcoming the above two difficulties, the experi- 
mental realization of our scenario amounts to controlling 
three parameters whose interplay underlies the physics 
of this work: the strength of an optical lattice s, the 
mean-field interaction between bosonic atoms griT f , arL d 
the strength of disorder R. All these quantities are ex- 
perimentally controllable using state-of-the-art technolo- 
gies. The interatomic interaction can be controlled in 
a very versatile manner via the technology of Feshbach 
resonances (40[ |. In the typical experiments to date, the 
values of ratio (JUtf/Eh range from 0.02 to 1 The 
depth of an optical lattice s can be changed from OEr to 
32Eji almost at will [4lJ . The above scenario for disorder 



can be realized in cold atomic systems using several meth- 
ods in a controlled way. They include applying optical 
potentials created by laser speckles or multi-chromatic 
lattices [42-45], introducing impurity atoms in the sam- 
ple [46| and manipulating the collision between atoms 
[47|] . Therefore, the phenomena discussed in this paper 
should be observable within the current experimental ca- 
pability. We emphasize here that the work presented in 
this paper is restricted to weak disorder and weak inter- 
atomic interaction. For further investigations in the pres- 
ence of stronger inter-atomic interaction or disorder, the 
pj^th-integral Monte Carlo simulation is a reliable method 

In summary, we explore the combined effects of weak 
disorder and a 2D optical lattice on the collective exci- 
tations of a harmonically trapped BEC. Accordingly, we 
extend the hydrodynamic equations for superfluid [8] for 
a 2D optical lattice to include the effects of weak disor- 
der. For small amplitude excitations, the corresponding 
equations are characterized by the disorder-depleted 
superfluid fluid density, a renormalized interatomic 
interaction and an effective mass capturing the effect of 
periodic potential. Using these hydrodynamic equations, 
we calculate the collective frequencies beyond the 
mean-field approximation. Our results reveal peculiar 
role of disorder, interplaying with optical lattice and 
atom-atom interactions, on the elementary excitations 
of a trapped BEC. In particular, we present detailed 
analysis for consequences of disorder on phonon propa- 
gation and lattice-modified surface modes. Conditions 
for possible experimental realization of our scenario are 
also proposed. 
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